Inelastic electron tunneling via molecular vibrations in single-molecule transistors 
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In single-molecule transistors, we observe inelastic cotunneling features that correspond energeti- 
cally to vibrational excitations of the molecule, as determined by Raman and infrared spectroscopy. 
This is a form of inelastic electron tunneling spectroscopy of single molecules, with the transistor 
geometry allowing in-situ tuning of the electronic states via a gate electrode. The vibrational fea- 
tures shift and change shape as the electronic levels are tuned near resonance, indicating significant 
modification of the vibrational states. When the molecule contains an unpaired electron, we also 
observe vibrational satellite features around the Kondo resonance. 

PACS numbers: 73.22.-f,73.23.-b,73.23.Hk 
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Electron tunneling is widespread throughout chem- 
istry and condensed matter physics. Electron transfer 
through molecules by nonresonant tunneling has long 
been known 0, 0], and tunneling electrons can interact 
inelastically with molecules, exciting vibrational modes. 
This is the basis of inelastic electron tunneling spec- 
troscopy (IETS)0,|3|, recently refined to probe groups 
of molecules via crossed wires 0, 0, and nanoporesQ, 
and single molecules via scanning tunneling microscopy 
(STM)^]. Similar inelastic, nonresonant tunneling oc- 
curs in single-electron transistors (SETs)0,^3, three- 
terminal devices in which a gate electrode allows in-situ 
adjustment of the energy levels. To date, this tuning has 
not been possible in the chemical systems examined by 
lETS. 

In this Letter we report inelastic cotunneling processes 
in single-molecule transistors (SMTs)|ll El El IH, 
and identify them with vibrational excitations of the 
molecules, as determined by Raman and infrared spec- 
troscopy. Tuning electronic levels near resonance reveals 
shifts in inelastic lineshapes and peak positions, suggest- 
ing significantly modified electron- vibrational coupling in 
this region. The vibrational features persist in the Kondo 
regime, indicating a complicated conduction process in- 
volving vibrational excitations of a many-body electronic 
system. 

A confined electronic system ("island") coupled by 
tunneling to source and drain electrodes is shown 
schematically in Fig. The energy to promote an 

electron from the highest occupied (HO) to the lowest 
unoccupied (LU) orbital is the single-particle level spac- 
ing, A. Ignoring spin, adding an electron to the island 
also requires additional Coulomb "charging energy", Ec, 
often approximated in SETs by the charging energy of 
a classical capacitor. In molecules, a manifold of vi- 
brational excitations is associated with each electronic 
state. In a three-terminal device, the island levels may 
be shifted via a gate potential. Fig. ^ maps the dif- 
ferential conductance, dlo/dYsD, of a generic SET as a 
function of source-drain bias, Vsd, and gate voltage, Vg, 
while Fig. QJ; shows possible electronic transport mecha- 
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FIG. 1: (a) Energetics of single-electron transistor, with A = 
single-particle level spacing, Ec = charging energy. Dashed 
levels indicate manifold of vibrational excitations in molecular 
devices, (b) Map o{ dlo/dVsD (brightness) of such an SET, 
vs. source-drain bias Vsd and gate voltage Va, showing one 
Coulomb transition from n to n+1 electrons on the island, (c) 
Conduction processes corresponding to the numbered points 
in (b): 1 = elastic cotunneling; 2 = inelastic cotunneling; 3 
= resonant conduction; 4 — resonant inelastic conduction; 5 
= Kondo resonant conduction. 



nisms. For low biases (Fig.^, 1,1') the average number 
of electrons on the island is fixed; transport is suppressed 
(Coulomb blockade) and can only occur by higher order 
tunneling through virtual states. An example of this in 
chemical electron transfer is "superexchange" 2|, and in 
SETs such processes are called "elastic cotunneling" jliJ] . 
At higher biases in the blockaded regime, "inelastic co- 
tunneling" via an excited virtual state (Fig.^:, 2) is pos- 
sible. For an excitation of energy E* , the opening of the 
inelastic channel results in a feature in d^Io/dVgjj at 
sVsD = E* . Inelastic cotunneling via electronically ex- 



2 



cited states has been seen in semiconductor 11] and car- 
bon nanotube SETslI^. Inelastic cotunneline; via vibra- 
tionally excited molecules is responsible for conventional 
IETST7], but has not been studied in three-terminal de- 
vices. lETS lineshapes are predicted to vary significantly 
depending on the energetics of the virtual states [l3|. 
and can be peaks, dips, or intermediate structures in 
d'^Io/dVgjj. At still higher source-drain biases (Fig.QJ;, 
3) Coulomb blockade is lifted leading to significant reso- 
nant conduction, while at still higher biases (Fig. 4) 
additional resonant conduction occurs when e V.qp is suf- 
ficient to leave the island in an electronicallv |l9l l2Ct I2H 
or vibrationallylH IH IH El excited state. 

When the effects of unpaired spins are included, the 
Kondo resonance becomes a possible conduction mecha- 
nism, and is detected as a sharp conductance peak near 
zero bias [23, 01 at temperatures low compared to a char- 
acteristic energy scale, fee Ik- The maximum peak con- 
ductance possible is 2e^//i, for a system with perfectly 
symmetric coupling between the island and the source 
and drain electrodes. Kondo physics has been observed 
recently in single-molecule transistors 13, 14, 15] . 

Single-molecule transistors open the possibility of 
examining inelastic cotunncling in individual, tunable 
chemical systems, in both the blockaded and Kondo 
regimes. We fabricate SMTs (Fig.l^J^) using an electromi- 
gration method that has been described extensivelvfl^ 
The source and drain electrodes are 15 nm 
Au films with 1 nm Ti adhesion layers, prepared by elec- 
tron beam lithography, e-beam evaporation, liftoff, and 
oxygen plasma cleaning. The gate oxide is 200 nm Si02, 
with a degenerately doped p+ Si (100) wafer as the un- 
derlying gate electrode. While this thickness of oxide 
guarantees relatively weak gate couplings due to simple 
geometric considerations, we routinely obtain excellent 
device reliability up to gate fields of 5 x lO^V/m with 
negligible leakage. 

The starting molecule, 1, (Fig. 13)) comprises a single 
transition metal ion (Co^"*") coordinated by conjugated 
ligands; the valence state of the ion may be controlled 
electrochemically (Fig. |2t). The structure of 1 has been 
verified by x-ray crystallography. Compound 1 under- 
goes loss of the (CN) moieties upon assembly on gold in 
tetrahydrofuran (THF), yielding the corresponding di- 
or tri-thiolate, which is bound covalently to surface Au 
atoms '255. Both 1 and its self-assembly have been char- 
acterized extensively by x-ray photoemission, ellipsome- 
try, and electrochemical methods [2^. 

A 2mM solution of 1 in THF is allowed to self-assemble 
on the electrode sets for 48 hours. The chip is then rinsed 
in THF, dried in a nitrogen stream, and placed in a vari- 
able temperature vacuum probe station (Desert Cryogen- 
ics) . This cryostat does not have magnetic field capabili- 
ties, but does permit characterization of many junctions 
at once. The junctions are partially broken to the kO 
level by electromigration at room temperature, cooled to 
liquid helium temperatures, and broken by further elec- 
tromigration into separate source and drain electrodes. 
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FIG. 2: (a) Diagram of typical single-molecule transistor fab- 
ricated by electromigration. (b) Structure of 1, the com- 
pound of interest, (c) Cyclic voltamrnograms showing re- 
versible change in Co ion charge state |2al . (d) Differential 
conductance of a transistor containing 1, at 5 K. White cor- 
responds to din /dVsD = 3 X 10~^ S. The zero-bias peak in 
the righthand charge state is characteristic of Kondo resonant 
conduction. 

At a variety of temperatures, we have measured Ib 
as a function of Vsd at Vg from -100 V to +100 V 
using a semiconductor parameter analyzer (HP4145B), 
with the source grounded and the drain electrode swept. 
We compute the differential conductance dlu / dVsD and 
d'^IoldYg^ as a function of Vsd and Vq by numerical 
differentiation. Spot comparisons between this a ppr oach 
and lock-in techniques show excellent agreement . To 
avoid artifacts wc identify features in d^In/dVgjj by 
comparing multiple data sets taken at the same or nearby 
gate voltages. 

We have examined 407 electrode pairs on 10 sepa- 
rate substrates, with statistics similar to those in pre- 
vious investigations '155 . Of the electrode pairs exam- 
ined we found 57 devices with no detectable current 
(electrodes too far apart for measurable conduction); 
166 with linear, nongateable current- voltage character- 
istics (likely no molecule present at the junction); 108 
with nonlinear but ungateable current-voltage charac- 
teristics (either molecules or metal nanoparticles with 
negligible gate coupling); and 76 nonlinear but signifi- 
cantly gateable current-voltage characteristics. From this 
last group, four were identified as single-electron devices 
based on unintentionally produced metal nanoparticles. 
This identification was based on the observation of many 
regularly spaced Coulomb blockade regions with typical 
electron addition energies less than 50 meV. In working 
devices, bias sweeps were restricted to \Vsd\ < 200 mV 
to minimize the chances of current-induced irreversible 
changes. Device stability is poor at high current den- 
sities and temperatures significantly above 4.2 K, likely 
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due to the atomic diffusion of the Au electrode material. 

There are 27 devices that cleanly display a single 
Coulomb degeneracy point, the vast majority of which 
show a resonance at zero Vsd in one charge state identi- 
fied as a Kondo resonance, as in Fig.|2i- Control devices 
with electrodes exposed to THF without 1 never display 
such conductance properties. As can be seen from the 
edges of the Coulomb blockade region, the electron ad- 
dition energy for the device in Fig. exceeds 100 meV; 
this is typical. The Kondo properties are very similar to 
those reported in earlier SMTs based on Co^+-containing 
complexes (iSj. Typical Kondo temperatures as inferred 
from the temperature dependence of the Kondo reso- 
nance height and the low temperature width of the reso- 
nance are ~40 K. The low temperature limit of the res- 
onance peak height is often reduced from the theoretical 
maximum value of 2e^//i, indicating asymmetric coupling 
of the molecule to the source and drain electrodes. A full 
discussion of Kondo physics in 1 and related molecules 
will be reported elsewhere. 

In 12 devices, the conductance in the classically block- 
aded region and/or outside the Kondo resonance is large 
enough to allow clean measurements of d^Io/dVgjy. In 
Figs. 13 we show maps of this quantity as a function 
of Vsd and Vg in two different devices at 5 K. We 
have indicated two prominent features within the block- 
aded (Kondo) regime with black arrows. Features in 
d^Io/dVgjj of opposite sign are symmetrically located 
around zero source-drain bias, consistent with features 
expected from inelastic tunneling. Some asymmetry in 
shape is unsurprising, given that the low peak conduc- 
tance in the Kondo regime for these devices indicates 
significantly asymmetric coupling of the molecule to the 
leads. 




FIG. 3: Maps of d'^lD/dVlo 

as a function of Vsd and Vg 
at 5 K for two devices. Smoothing window in Vsd is 5 mV. 
Brightness scales are -8 x 10"^ A/V (black) to 3 x lO"'^ A/V 
(white), and -2 x 10"^ A/V (black) to 2 x lO"'^ A/V, re- 
spectively. The zero-bias features correspond to Kondo peaks 
in dlo/dVsD- Prominent inelastic features are indicated by 
black arrows. In both devices, when the inelastic features ap- 
proach the boundaries of the Coulomb blockade region, these 
levels shift and alter lineshape (white arrows). Black dashed 
line in left map traces an inelastic feature across across the 
boundary and into the Kondo regime. 

The d'^Io/dVgjj features in the blockaded region occur 
at essentially constant values of Vsd until Vq is varied 
such that the feature approaches the edge of the block- 
aded region. This constancy in Vg has been identified 



as a feature of inelastic cotunneling in semiconductor'll^ 
and nanotube^^ single-electron devices. The changes in 
these features as the boundaries of the Coulomb stabil- 
ity region are crossed indicate that the inelastic processes 
are native to the SMT itself, and not due to some parallel 
conduction channel. The inelastic modes occur at ener- 
gies low compared to A (~100s of meV), implying that 
the modes being excited are unlikely to be electronic. By 
examining the d^In/dVgjy data for features at fixed Vsd 
as a function of Vg, we identified a total of 43 candi- 
date inelastic features at Vsd < 120 meV, ranging in 
apparent width from 5 mV to 20 mV at ~5 K. The nar- 
rowest inelastic features clearly show a broadening as the 
temperature is elevated above 20 K. Further studies will 
seek to compare this broadening with standard theoreti- 
cal treatments of linewidths in lETS. 

Fig. 01 shows a histogram (bin = 1 mV) of the Vsd 
positions of those features, from detailed examination of 
O^Id I dVgjy vs. Vsd data. The lower panel shows Raman 
(Stokes) peak positions (data taken on a powder of 1) and 
IR absorption peak positions of 1 (in pellet form, blended 
with KBr) j2^] . The correlations between the spectra and 
inelastic features in d^lD/dVgjj confirm that vibrational 
inelastic cotunneling processes arc at work in the molecu- 
lar transistors. This is consistent with the observation of 
lETS in single molecules by STM^g. The fact that differ- 
ent SMTs exhibit different subsets of features is a natural 
consequence of the sensitivity of lETS to the nanoscale 
structure of the molecular junction 27'|. Comparisons 
with Raman of 1 and related compounds (e.g. ligands 
without -SCN; hgands without Co^+) suggest tentative 
peak assignments. The peaks near 36 meV and 55 meV 
are likely the Co-N stretchj^^, while that near 44 meV 
is tentatively a -SCN deformation[2^ (implying that not 
every S is bonded to an Au surface atom, consistent 
with x-ray photoemission dat al2al ). The 24 meV peak 
is tentatively the Au-S bond|29|: this claim is further 
supported by the appearance of a ^28 meV Raman peak 
in films of 1 assembled on Au electrodes, while such a 
peak is absent in Raman spectra of both unassembled 1 
and bare electrodes. 

As indicated by white arrows in Fig. |31 inelastic 
features shift and change near the boundary of the 
Coulomb blockade region. These changes as a function 
of Vg are qualitatively consistent with recent theoret- 
ical expectations 10, though detailed modelling would 
be required for quantitative comparisons. The energy of 
the virtual cotunneling states is shifted by Vg, altering 
the (complex) amplitude for that process. Changing the 
cotunneling amplitude compared to quasi-elastic (with 
the virtual excitation and readsorption of a vibrational 
quantum) or direct source-drain tunneling is predicted to 
severely alter lineshapes. These effects cannot be exam- 
ined in two-terminal tunneling structures, which lack the 
ability to shift the relevant virtual states. We note that 
occasionally inelastic features shift significantly in energy 
as the blockade boundary is approached (e.g. the 24 meV 
features in Fig.|2Ileft)). This would be consistent with a 
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FIG. 4: (top) Histogram of \Vsd\ positions of 43 features in 
dI%/dV§D that are persistent at constant Vsd over at least 
a 10 V gate voltage range, for the 12 samples discussed. Bin 
size is 1 meV. Width of actual features is at least 5 meV, with 
position indicating feature center, (bottom) Raman (black) 
and IR (gray) peak positions taken of 1 in the solid state at 
room temperature. The indicated Raman peak near 28 meV 
is only present in films of 1 self-assembled on Au electrodes. 

modified (dressed) electron- vibrational coupling near the 
electronic resonance. 

We observe that the lETS features persist into the 
Kondo regime in the form of satellite features paralleling 



the zero-bias Kondo resonance. Often specific features 
may be traced from the blockaded region into the Kondo 
region, an example of which is indicated by the dashed 
line in Fig. I^fleft). Through the transition the inelastic 
feature lineshapes change, as do their intensities. Such 
inelastic satellites paralleling the Kondo resonance have 
been suggested Dreviouslv |l4i . and the consistency in 
feature position in our data lend strong support to this 
idea. Recent experiments in semiconductor devices (SOj 
have demonstrated the existence of satellite Kondo peaks 
when the Kondo system can interact inelastically with 
photons of a well defined energy. In the molecular tran- 
sistor case, the inelastic exchange occurs with the vibra- 
tional quanta of the molecule, demonstrating a quantum 
mechanically coherent coupling between the electronic 
many-body Kondo state, and the mechanical resonances 
of the molecule. These Kondo measurements illustrate 
how the tunability of single-molecule transistors permits 
the examination of vibrational processes as a function of 
the energetics of the electronic levels, a study not possible 
in standard two-terminal devices. 
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